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locus corresponding to coil-globule transition clearly approaches the real axis as the chain 

length increases, and the transition temperature could be estimated by finite-size scaling. The 

form of the logarithmic correction to the scaling of the partition function zeros could also be 

obtained. The other locus does not show clear scaling behavior, but a supplementary analysis 

of the specific heat reveals a first-order-like pseudo-transition. 
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Introduction 



Conformational transitions of a polymer have been a subject of interest for many years. 1-32 The 
coil-globule transition of a homopolymer, also called the collapse transition, has been most widely 
studied. A polymer chain in a dilute solution is influenced by both hydrophobic interactions be- 
tween the monomers and the excluded volume effect, which cancel each other at a special tem- 
perature T = Tq. The polymer chain adopts expanded and compact conformations for T > Tq 
and T < Tq respectively, with T = Tq being the temperature of coil-globule collapse transition. 
T > Tq region of the polymer system can be mapped into a magnetic system where the number 
of components is formally set to zero, with the infinite chain length corresponding to the critical 
temperature^. 4 - 6 Then T = Tq is a tricritical point where both of the quadratic and the quartic terms 
in the corresponding field theory vanish. 6 

The radius of gyration (or the end-to-end distance) R of a polymer chain with N monomers, 
near T = Tq, is generally expressed by the scaling theory, 6 

(R 2 )~N 2v <f(TN*), (1) 

where the exponent v t represents the geometrical properties of a polymer at the tricritical point, and 
the exponent 0, called the crossover exponent, 33,34 measures how rapidly the system undergoes the 
transition as the temperature T approaches the tricritical temperature Tq . The reduced temperature 
is defined as T = \T — Tq \ /Tq, scaling as x ~ as T — > Tq. The scaling function f(x) behaves 
as follows— 

' x (6/(d+2)-2 Vt )/<t, if^oo 



/(*) 



const. if x — > 0, (2) 



'The temperature of the magnetic system is conjugate of the polymer chain length Af, and it has nothing to do with 
the temperature T of the polymer model. 
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Since the upper tricritical dimension is three, the tricritical exponent can be obtained from the 
mean field theory as § = 1/2 and v t = 1/2, but the scaling form (OQ) is expected to have logarithmic 
corrections. 5 

In addition to the coil-globule collapse transition which is of second order in three dimensions, 3 
a first-order transition from liquid-like globule to a solid-like phase is expected at a temperature 
lower than Tq,—~— which has been the subject of active research recently but less well understood 
than the coil-globule collapse transition. 

In this work, we enumerate the number of all possible conformations on a simple-cubic lattice 
up to chain length 24. In particular, we study the zeros of the exact partition function, which are 
much more sensitive indicators of phase transition than the real- valued quantities such as specific 
heat. We observe that for chain length longer than 16 there are two distinct loci of zeros, suggesting 
the existence of both the coil-globule collapse transition and liquid-solid melting-freezing transi- 
tion. The exact partition function zeros of three-dimensional lattice polymer have been studied for 
simple 12 and face-centered 8 cubic lattices for chain lengths up to 13 and 10, respectively, but these 
lengths were too short to reveal two distinct loci, and the results for longer chains up to 31 12 were 
based on Monte Carlo sampling with limited accuracy. 

Although the chains whose conformations are exhaustively enumerated in this work are still 
much shorter than those studied with Monte Carlo samplings, the exactness of the calculation 
enables us to extrapolate the finite N data to obtain the result in the limit of infinite chain length 
with a reasonable accuracy. The chains were still too short to extract the large N behavior of the 
inner locus corresponding to the melting-freezing transition, but the tricritical collapse temperature 
Tq could be obtained from the scaling behavior of the outer locus. The form of the logarithmic 
correction to the scaling of the partition function zeros could also be obtained. Additional analysis 
was performed on the exact specific heat, which shows the first-order-like pseudo-transition near 
the inner locus, in accordance with the previous result based on a chain-growth sampling. 26 
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The Number of Conformations 

Conformations of a polymer chain with N monomers are modelled by self-avoiding walks with 
N— 1 steps on a simple-cubic lattice. The position of a monomer i is expressed as r* = (a, b, c) with 
integer values of the coordinates a, b, and c. The coordinates satisfy the constraints |r, — r, + i | = 1 
and r ; 7^ r ; for i 7^ j due to the chain connectivity and the excluded volume effect. 
We consider the Hamiltonian with the nearest-neighbor interaction: 



and £ is set to a positive value to incorporate the attractive interaction between the monomers. 

Assuming that the polymer chain has an intrinsic direction, the conformations with reverse 
labels i ■<->■ N — i + 1 for all (i = 1, 2, • • • ,N) are considered distinct. For a generic conformation, 
the rigid rotations and reflections form an 48-fold symmetries in three dimensions. Exceptions 
are the cases of lower-dimensional conformations embedded in higher dimensional spaces. In 
three dimensions, 24-fold and 6-fold symmetries exist for the planar and linear conformations, 
since they are invariant under transformation perpendicular to the underlying plane and straight 
line. Therefore reduced numbers of conformations, where conformations related by symmetry 
are counted only once, are computed in order to prevent the waste of computational resources.— 
Since the energy © depends only on the number of inter-monomer contacts K, we classify the 
conformations according to the value of K. The number of conformations Q.^ (K) for a given 
contact K in d dimensions, with discrete rotations and reflections considered distinct, is obtained 



^ = -e£A(r ; ,r 7 ), 



(3) 



where 




(4) 
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from the reduced number of conformations (0^ (K) as — 

nW(K) = 2coW(K) = 28 Kfi , 

a^(K) = S(0^(K)+4(0^(K), (5) 
Qf?) (K) = 48G)( 3 ) + 24w( 2 ) (K) + 6(0^ (K). 

From here on, we drop the dimension index and use Q(K) to denote Cl@)(K). The number of 
conformations Q,(K) for Af = 24 is presented in Tabled] 

Partition Function Zeros in the Complex Temperature Plane 

Partition function zeros have been the subject of interest as a sensitive indicator of a phase transi- 
tion.^ - — i 3 - 7 -"^ 3 - Partition function zeros were introduced by Yang and Lee in the complex fugacity 
plane of a fluid system and the complex magnetic-field plane of the nearest-neighbor Ising fer- 
romagnet (Yang-Lee zeros), to study the phase transition driven by the fugacity or the magnetic 
field. 37 Later, Fisher 38 used the partition function zeros in the complex temperature plane (Fisher 
zeros) of the square-lattice Ising model, to study the temperature driven transition. In the ther- 
modynamic limit, the locus of zeros forms a continuous curve which crosses the real axis if a 
transition exists. Thus, the theory of partition function zeros provides the explanation on how the 
partition function, which is an analytic function of thermodynamic parameters at a finite size, ac- 
quires the singularities necessary for a phase transition in the thermodynamic limit. In the case 
of Fisher zeros, the transition temperature in the thermodynamic limit is the intersection point of 
the locus of zeros with the real temperature axis. Therefore, the conjugate pair of zeros closest to 
the positive real axis, called the first zeros, determine the leading singular behavior of the partition 
function. Since the behavior of the first zeros can be analyzed separately from the other zeros, 
the phase transition can be analyzed more accurately by computing the partition function zeros 
than studying real- valued quantities such as the specific heat which includes the effect from all the 
zeros. 
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In the case of the lattice polymer of our interest, the partition function is written in terms of 
Q{K) as 

Z(y) = ^e-^ = ^a(K)y- K , (6) 

{CT} K 

where {0} denotes a sum over all possible conformations, /3 = l/kgT, and y = exp(— /3e) which 
ranges between and 1 for positive T; y = when T = and y = 1 when T — > °°. Since the 
maximum number of contacts Km is finite for a given chain length N, y KM Z{y) is a polynomial of 
order Km in y. Therefore the partition function zeros in the complex y plane can be obtained by 
solving the polynomial equation y KM Z(y) = with MATHEMATICA. 

Figured] shows the partition function zeros in the complex temperature plane for several values 
of chain lengths. The result for N = 13 (Fig. [Ha)) agrees with the previous result obtained by 
Finsy et al. 12 by exact enumeration after a change of variable, but the chain length is too short to 
reveal two distinct loci of zeros. The partition function zeros for 14 < N < 31 were studied in the 
same reference by Monte Carlo sampling, but again two distinct loci were not visible. One of the 
reasons is that the zeros in the complex temperature plane of z = e^ £ were investigated there instead 
of y = e~$ £ = 1 jz as in the current work. In the former case, the temperature range of interest T > 
corresponds to the real line with z > 1, with T = corresponding to z = +°°. Therefore, features 
relevant to low temperature behavior of the system may be easily missed if only a finite region near 
origin is considered. In contrast, T > corresponds to the segment < y < 1 in the y plane, and 
consequently all of the physical regions can be examined with ease. 

The exact zeros plotted in the y plane clearly exhibit the splitting of the locus into two distinct 
branches already for N > 16 (Fig.[TJb) and (c)). Since there are two visually distinct loci of zeros, 
we will select the first zeros for each locus and analyze them separately in the following sections. 

The Outer Locus and the Coil-Globule Collapse Transition 

The first zeros of the outer locus approach the real axis suggesting a non-trivial transition in the 
thermodynamic limit (Fig. |2ta) (open circles)). From the scaling of the partition function whose 
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form is similar to (0Q>, the first zeros y\ are expected to scale as:— >22 

yi(N)-y e ~N-*, (7) 

where ye = exp(— e/ksTe). In three dimensions corresponding to the upper tricritical dimension, 
the scaling form CQ) is modified by logarithmic corrections,^*^ and consequently © is expected 
to be modified to be of the form 

yi(N)-y e ~N-*(\ogN)- x . (8) 
The power of logarithmic correction A can be estimated from the imaginary part of d8}: 

lm[ yi (N)] ~Ar4>(l O gA0~ A , (9) 

which is rewritten as 

Jm\yi(N)]N* ~ (logA0~ A . (10) 

Taking the logarithm, we get 

log(Im[yi ~ - A loglogN + constant. (11) 

In order to remove the constant term, the difference of ([Til for neighboring chain lengths is taken: 

Am[vi(iV+2)](Ar+2)n „ / log(* + 2) \ 

10 H ImbK^V)]^ Al0g { logN )> U2) 

where we used N + 2 instead of N + 1 because even and odd number of chain lengths are expected 
to exhibit distinct scaling behaviors.— From (fT2)) the finite-size estimation of A is obtained as: 

log(Im[yi(jy + 2)](jy + 2)»)/Im[yi(A0]A^) 
1 ] ~ log(log(jV + 2)/logA0) 1 ^ 
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with = 1/2. We applied Bulirsch-Stoer (BST) method 44 for the data of even N with 12 < /V < 24 
to estimate the value in the limit N — > °°. We obtained 

A =0.642(28), (14) 

where the error is estimated by examining the robustness of the extrapolated value with respect to 
perturbations of the data points.— The data points are chosen to maximize this robustness. The 
estimated value (fT4l) exhibits remarkable agreement with 

^ — tj =0.636346, (15) 

appearing in the scaling of the Boyle temperature 7g where the second virial coefficient van- 
ishes:^^ 

T B (N)-T e ~Ar 1/2 (logA0~ 7/U . (16) 

We then estimate the collapse temperature Tq by taking the real parts of the scaling relations © or 
®. Applying the BST extrapolation for the data of even N with 14 < N < 24, the collapse tem- 
perature is obtained as ye = 0.7185(94) (fcgTg/e = 3.03(12)) without the logarithmic correction, 
and y e = 0.7653(174) (k B T e /e = 3.76(32)) with the correction factor of (log A^)" 7 / 11 (Fig. Efo)). 
As shown in Table [2l the agreement with previous results is better when the logarithmic correction 
is included. It is to be noted that although the maximum length of chain 24 considered here is 
much less than those in Monte Carlo studies which is as long as 32000, 26 the exactness of our data 
enables us to perform extrapolation to large /V with reasonable accuracy. 

We also note that the intersection of the locus of zeros with the real axis is estimated to be 
at a temperature higher than ksTg/e = 1.81(2) obtained by fitting the Monte Carlo data of par- 
tition function zeros to a polynomial curve.— It was conjectured in the same reference that the 
corresponding transition is of first-order, which was later conjectured to be a low temperature 
melting-freezing transition. 3 However, we found no evidence that the transition is of the first or- 
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der. Since there is a separate inner locus corresponding to a transition at a lower temperature, and 
the intersection of the outer locus with the real axis in the thermodynamic limit is now estimated 
to be at 3.76 rather than 1.81, the outer locus seems to be the one corresponding to the coil-globule 
collapse transition, which is well established to be a second-order transition.- - - 



The Inner Locus and the Melting-Freezing Transition 

The scaling of inner locus is much worse than that of the outer locus in this range of chain lengths. 
The distribution of the first zeros for various chain lengths is rather irregular (Fig. 0a)), and the real 
part shows oscillatory behavior as a function of TV -1 / 2 (Fig.|2£b)). Therefore no reasonable estimate 
of the transition temperature in the thermodynamics limit could be made. This irregular behavior 
of the melting-freezing transition point is also in agreement with the previous result obtained by 
computing the specific heat with chain-growth sampling for /V up to 125. 26 This irregularity is 
understandable, as explained in the same reference. The effect of melting-freezing transition which 
is expected to survive in the limit of large /V limit, and that of the excitation pseudo-transition 
which appears only near a special value of N called magic numbers, are intermixed to give a rather 
complex behavior. 

To complement the result from the partition function zeros, we performed additional analyses 
by computing the exact specific heat per monomer, 

Cv(T,N) 1 dE /3 2 <9 2 lnZ (\ny) 2 , 2 2 

-NkT = M- B dT = -N^ = — ( {K) - {K) )' (17) 

which is plotted in Fig. Oa) as a function of y for several values of N. The effect of the melting- 
freezing and excitation transitions are mixed to manifest themselves as one prominent peak. In 
accordance with the previous result from a chain growth sampling, 26 we find that due to the effect 
of the excitation transition at the magic numbers N c = 8, 12, 18, ■ ■ ■ satisfying N c = L? or N c = 
L 2 (L± 1), the peak becomes sharper as /V is increased from N c — 1 to N c and becomes flat again at 
N = N c + 1 (Fig.[3ta)). As to be expected, the positions and the values of the peak exhibit irregular 



behavior as functions of the chain length (Fig. |3{b)), similar to the result from the chain growth 
sampling. 26 The peak positions are in the range 0.10 < y pea k < 0.25, much lower than the coil- 
globule collapse transition temperature y c = 0.7653(174) obtained from the analysis of the outer 
locus of the partition function zeros, but much closer to the inner locus of zeros, whose real values 
are distributed mostly in the range 0.0 < y < 0. 1 (Fig. Ufa)). 

In the microcanonical formalism, the first-order phase transition is signaled by the existence 
of a region of energy with ^§ > where S(E) = log^E 1 ) 2 , which is invisible in the canonical 
formalism.^^ 6 - For the system with discrete values of energy, this condition translates into the 
condition that there exists i with 



S(E i+1 ) -2S(E i )+S(E i ^ l )>0, (18) 

where i labels the energy values in the ascending order. These conditions are equivalent to the 
existence of j8o near the transition temperature such that 

e- hE a(E) (19) 



has at least two peaks corresponding to distinct phases. We see that the chains with N > 16 as well 
as even value of N > 8 satisfy these properties. We plot the two peaks of e^^ )E Q.(E) in Fig. @]for 
the largest magic number (N = 18) and the largest number (N = 24) among the chain lengths we 
studied, where the values of /3o were adjusted for each chain so that the heights of the peaks are 
the same. These features clearly shows the first-order-like nature of the pseudo-transition, and are 
also consistent with similar results from the chain-growth sampling. 26 

In contrast to melting-freezing and excitation transition, the effect of the coil-globule collapse 
transition is not readily visible in the specific heat (Fig. |3]). This is due to the fact that the collapse 
transition is of the second order, where the specific heat is continuous and only its derivative is 
discontinuous or divergent. 

2 This corresponds to a negative value of microcanonical specific heat C m i cro = — j3 2 ( J-^gT^ J < 0- 
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Discussions 



We studied conformational transitions of a polymer by exhaustively enumerating the number of all 
possible self-avoiding walks on a simple-cubic lattice up to chain length 24. Although the lengths 
of the chain are much smaller than those studied by Monte Carlo samplings, the strength of our 
result is that they are obtained from exact enumeration and hence contain no errors, enabling us to 
make extrapolation to infinite chain length with a reasonable accuracy. 

Furthermore, by studying partition function zeros in the complex temperature plane, we could 
obtain information which is not readily available in the real-valued quantities such as specific 
heat. We observed two distinct loci of partition function zeros in the complex temperature plane, 
suggesting the existence of coil-globule collapse transition and liquid-solid melting-freezing tran- 
sition. 

From the finite-size scaling of the first zeros of the outer locus with mean field crossover expo- 
nent = 1/2, the scaling form of the first zero with logarithmic correction factor was conjectured 
to be y — yg ~ /V~ 1//2 (log/V)~ 7//11 as in the case of the scaling of the Boyle temperature. The col- 
lapse transition temperature was estimated to be 1c b Tq /e = 3.03(12) and KbTq /e = 3.76(32) with 
and without logarithmic correction, respectively. The result shows better agreement with previous 
results in the presence of the logarithmic correction, suggesting that the conjectured form of the 
logarithmic correction to the scaling of the first zeros is indeed correct. 

The results for the collapse transition indicate an additional advantage of studying the partition 
function zeros. Being a second-order transition where the specific heat is finite and continuous, 
the signal for the transition cannot be easily detected by examining the peak of the specific heat 
(Fig. Sa)). On the other hand, since the partition function zeros is due to the singularities of the 
partition function regardless of the order of the transition, the signal for the transition is clearly 
visible as a locus of zeros approaching the real axis as the chain length increases (Fig. \T\ and |2]), 
making them an indispensable tool for studying phase transitions. 

The behavior of the inner locus was not regular enough for the current chain lengths to esti- 
mate the melting-freezing temperature in the limit of infinite chain lengths, due to mixture of the 
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effect from the finite-size excitation pseudo-transition. The existence of the energy region with 
the negative value of microcanonical specific heat, which manifests itself as double peaks in the 
density of states multiplied by Boltzmann factor near transition temperature, shows a clear sign of 
first- order-like nature of the pseudo-transition. 

The excitation transition is analogous to the folding transition of HP protein, which appears 
only for a particular sequence. 47 ^2 It would be interesting to compare the transition behaviors of 
designable proteins, random heteropolymers and homopolymers, all on the same lattice. 

It is to be noted that although partition function zeros can be used to investigate the existence 
and the property of the phase transition, the information on the nature of the phases themselves 
cannot be obtained. Various geometric parameters of conformations contributing the phases of 
interest should be analyzed in order to confirm that the transitions we observe are indeed the coil- 
globule collapse and solid-liquid melting-freezing transitions. 
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Table 1 : The number of conformations Q.(K) as a function of the number of contacts K for N 



K 


Q(K) 





238306751550942 


1 


601441550088000 


2 


856234452257592 


3 


919771036344192 


4 


821203501326936 


5 


642255091228800 


6 


455089815998760 


7 


298905402843360 


8 


184343422767744 


9 


107615281912368 


10 


59739246931968 


11 


31649589839232 


12 


16004806431576 


13 


7677745597008 


14 


3470790178464 


15 


1454509923624 


16 


559820945808 


17 


190765562640 


18 


57066241104 


19 


13933700784 


20 


3113368896 


21 


477160080 


22 


30437280 


23 


12554256 


Total 


5245988215191414 
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Table 2: The values of the coil-globule collapse temperature ksTQ/e of a lattice polymer on a 
simple-cubic lattice with nearest neighbor interaction, obtained in the current work (first two lines), 
are compared with previous results on the same model. 



ivieinou 


AI 

ly max 


K-bIq/ t 


Exact FrZ (without (log/VJ ' correction) 


24 


3.0.3(12) 


Exact PFZ (with (log/vj ' correction) 


24 


3.76(32) 


RCG^- 


2000 


3.64 ~ 4.13 


MC Ft 1 / 


31 


1 O 1 /1\ 

1.81(2) 


MC 13 


40 


1 "~7 1 1\ 

3.714(H) 


* /r/~i 14 

MC 


299 


4 


MC 


1024 


3.713(7) 


SS— 


250 


3.65(8) 


MC 


100 


2.972(6) 


MC— 


5000 


3.721(6) 


CBV 19 


200 


3.45 


WCM 


1200 


3.598(54) 


PERM 21 


10000 


3.724 


PERM 22 


2048 


3.717(2) 


MC^ 3 - 


1000 


3.71(1) 


BGY 24 


600 


3.745 


MC^ 


16000 


3.71(1) 


Improved PERM 26 


32000 


3.72(1) 



2 PFZ: partition function zeros 

2 RCG: Rosenbluth chain growth 

2 MC: Monte Carlo 

2 SS: scanning simulation 

2 CBV: configurational-bias-vaporization 

2 PERM: pruned-enriched Rosenbluth method 

2 BGY: Born-Green- Yvon integral equation 
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Figure 1: Distribution of the partition function zeros in the complex temperature (y = e $ £ ) plane 
for (a) N = 13, (b) N = 16, and (c) N = 24. 
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Figure 2: (a) The first zeros of outer locus in Fig. [T]for /V = 10, 11, 12, ■ • ■ , 24 are plotted as open 
circles, and those of inner locus for N =16, 17, 18, ■ ■ • , 24 as open squares. Since the distribution 
of zeros is symmetric with respect to the real axis, only the first quadrant of the complex plane is 
shown, (b) Values of the real part of the first zeros of the outer locus are shown as a function of 
/V~ 1//2 (solid circles) and /V~ 1//2 (ln/V)~ 7//11 (solid triangles) for even N with 14 < N < 24. From 
the extrapolation to N — > °° = 0), the collapse temperature is obtained as ye = 0.7185(94) 
without a logarithmic correction (the open circle with an error bar), and y$ = 0.7653(173) with 
the correction factor of N~ l / 2 (lnN)~ 7 / 11 (the open triangle with an error bar). On the other hand, 
those of the inner locus forN= 16, 17, 18, • • • , 24 (open squares, plotted as the function of N~ 1 / 2 ) 
show irregular behavior and cannot be extrapolated. 
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Figure 3: (a) Curves of the specific heat per monomer as a function of y, for N near the magic 
number 18 = 2- 3 2 . The peak rises sharply as N is increased from 17 to 18, but becomes flat 
again as it is increased to 19, due to effect of the excitation pseudo-transition at /V = 18 (b) Peak 
temperatures of the specific heat for 8 < N < 24 as a function of chain length N, showing an 
irregular pattern within 0. 10 < y pea k < 0.25. 
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E (N=24) 




E(N=18) 

Figure 4: The density of states multiplied by Boltzmann factor e~P ()E Q.(E) as a function of E, with 
/3o near the transition temperature, forN= 18 and 24. The values of /3o are adjusted for each N 
to make the heights of the two peaks equal. The double peaks show a clear sign of first-order-like 
transition. 
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